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Abstract.  A  mechanistic  model  that  interprets  the  transition  in  oxidation  behavior  of  zirconium 
diboride  as  the  temperature  is  varied  from  600°C  to  2500°C  is  presented.  Available  thermodynamic 
data  and  literature  data  for  vapor  pressures,  oxygen  permeability  in  boria,  and  viscosity  of  boria 
were  used  to  evaluate  the  model.  Three  regimes  and  the  temperatures  of  transition  between  them 
were  identified.  In  the  intermediate  temperature  regime,  viz.,  1000°C  to  1800°C,  good 
correspondence  was  obtained  between  theory  and  experiments  for  weight  gain,  recession,  and  scale 
thickness  as  functions  of  temperature  and  oxygen  partial  pressure.  In  this  regime,  the  rate-limiting 
step  is  the  diffusion  of  dissolved  oxygen  through  a  film  of  liquid  boria  in  capillaries  at  the  base  of 
the  oxidation  product.  At  lower  temperatures,  an  external  boria  scale  forms,  but  it  was  not  found  to 
contribute  significantly  to  oxidation  resistance.  Comparison  with  literature  data  on  recession  is  very 
good,  but  weight  gain  is  predicted  to  be  higher  than  experimentally  observed  unless  flow  of  viscous 
boria  is  included.  At  higher  temperatures,  the  boria  is  lost  by  evaporation,  and  the  oxidation  rate  is 
limited  by  diffusion  of  molecular  oxygen  through  the  capillaries  between  nearly  columnar  blocks  of 
the  oxide  MO2.;  this  regime  is  soon  followed  by  a  rapid  acceleration  of  recession  due  to 
vaporization  of  the  oxide  MO2  itself. 


Introduction 


Diborides  of  refractory  elements,  either  in  pure  or  in  composite  fonn  with  SiC  additions,  are  known 
to  be  resistant  to  oxidation  and  erosion  at  high  temperatures  [1-8].  More  recent  work  has  been 
prompted  by  the  possible  application  in  hypersonics[9-17].  However,  little  or  no  work  has  been 
reported  on  quantitative  models  of  oxidation  of  this  class  of  materials.  A  few  recent  works  focused 
on  thermodynamic  studies[15,  18]  which  have  identified  the  key  solid  and  gaseous  oxidation 
products.  Modeling  the  kinetics  of  oxidation  of  this  class  of  materials  has  been  the  focus  of  our 
research  program. 

In  the  first  study,  we  reported  on  a  model  for  the  oxidation  of  these  materials  in  the  intermediate 
temperature  regime,  viz.,  ~ I OOOC  to  1800C[19].  In  this  report,  we  examine  the  mechanisms  that 
operate  at  temperatures  below  1000C,  and  at  temperatures  above  1800C.  We  have  retained  the 
essential  elements  of  our  prior  model,  but  extend  it  to  include  effects  that  characterize  the  transition 
to  the  low-temperature  and  high-temperature  regimes. 

1 
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We  start  with  a  brief  overview  of  the  model  for  the  intermediate  temperature  regime.  Then 
equations  are  developed  that  describe  the  low  and  high-temperature  regimes,  before  comparing  the 
model  to  experimental  data.  The  experimental  data  on  these  materials  come  from  a  variety  of 
sources,  whose  processing  methods  have  varied  with  technological  advances  over  the  past  several 
decades.  The  materials  have  also  been  studied  under  a  variety  of  different  environmental  conditions, 
temperature  regimes  and  using  different  methods.  Thus  the  model  has  to  be  robust  to  include  the 
effects  of  all  experimental  variables  to  allow  for  a  fair  comparison  with  literature  data. 


Model 

The  key  elements  of  the  model  are  described  in  detail  in  [19].  Here  we  state,  in  brief,  the  key 
assumptions  of  the  model.  Fig.  1  shows  a  schematic  of  the  assumed  oxidation  products  and  their 
morphology  at  steady  state.  The  assumptions  are  based  on  reported  experimental  observations  as 
noted  in  [19].  At  low  temperatures,  i.e.,  below  about  1000°C,  the  scale  is  made  of  columnar 
zirconia  under  an  external  layer  of  liquid  boria.  The  boria  wets  the  zirconia  and  is  continuous  across 
the  zirconia  scale.  Oxygen  diffuses  through  the  external  boria  layer  and  then  through  the  tortuous 
pathway  between  zirconia  grains  to  reach  the  substrate.  The  boria  and  zirconia  that  form  at  the 
ZrEF/ZrCF  interface  continually  increase  the  zirconia  layer  thickness  and  the  external  boria 
thickness.  The  zirconia  that  forms  is  stable  and  is  resistant  to  oxygen  penetration  due  to  very  low 
ambipolar  diffusion  kinetics  (see  appendix  of  [19]).  However,  the  boria  either  flows  (viscous)  or 
evaporates  at  the  free  surface  (boundary  layer  diffusion  limited)  tending  to  reduce  the  external  boria 
thickness.  At  some  temperature,  around  1000°C  (depending  on  the  actual  dynamic  gas  flow 
conditions),  the  boria  evaporates  so  fast  that  the  scale  can  no  longer  keep  an  external  layer  of  boria. 
At  this  temperature,  the  oxidation  transitions  to  an  intermediate  temperature  regime,  where  the  boria 
is  still  present  in  the  scale,  but  only  within  the  zirconia  pores,  which  forms  a  continuous  pathway. 
Finally,  diffusion  in  pores  is  assumed  to  follow  Fick's  law,  and  kinetic  laws  are  derived  for  planar 
symmetry. 


Very  High  Temp  >  ~ 1800  C 


Fig.  1 .  Schematic  sketch  of  the  oxidation  products  and  morphology  assumed  in  the  model. 
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At  intermediate  temperatures,  the  porous  channels  within  zirconia  are  partially  filled  with  boria; 
thus  gaseous  oxygen  must  now  diffuse  through  part  of  the  zirconia  pore  channel,  and  then  penneate 
through  the  boria  in  the  pore  channel  the  rest  of  the  way  to  reach  the  substrate.  The  model  assumes 
that  a  quasi-steady  state  is  reached  where  the  amount  of  boria  formed  is  balanced  by  the  amount  of 
boria  that  evaporates,  which  detennines  the  depth  to  which  the  zirconia  pores  are  filled  with  boria. 

In  this  temperature  regime,  the  ambient  condition  is  taken  to  be  a  perfect  sink  for  boria.  The 
equations  governing  this  set  of  mechanisms  give  a  constant  ratio  of  the  boria  thickness  to  zirconia 
thickness  at  any  given  temperature[19].  As  the  temperature  increases,  this  ratio  decreases  until  a 
temperature  is  reached  where  the  boria  evaporates  and  diffuses  through  the  porous  channel  as  fast 
as  it  is  formed.  Thus  the  mechanism  evolves  towards  a  high  temperature  regime.  In  the  high 
temperature  regime,  there  is  no  liquid  boria;  the  scale  is  a  dry  porous  zirconia.  The  oxidation  is 
limited  by  the  rate  of  oxygen  diffusion  through  the  porous  channel.  In  this  temperature  regime, 
ZrCF  also  evaporates  significantly,  which  enhances  recession.  The  evaporation  kinetics  will  depend 
on  the  actual  dynamic  flow  of  gases. 


Model  Equations 

The  equations  describing  the  intermediate  temperature  regime  can  be  found  in  Ref. [19].  Equations 
that  relate  to  mass  balance,  gaseous  diffusion,  oxygen  permeation  rates,  and  thermodynamic  data 
are  the  same  for  the  low  temperature  and  high  temperature  regimes  and  hence  will  not  be  repeated 
here;  the  reader  is  referred  to  ref.  [19]. 

Low  temperature  regime.  A  schematic  sketch  of  the  model  for  low  temperature  oxidation  kinetics 
is  shown  in  Fig.  2.  The  length  of  the  zirconia  layer  is  L,  and  the  thickness  of  the  external  boria  is 
hext.  The  oxygen  partial  pressures  are  denoted  by  Po2-  We  use  superscripts  to  specify  the  value  of 
parameters  at  the  interfaces.  The  superscripts  s,  zb  and  a  refer  to  the  substrate/scale  interface,  the 
zirconia/boria  interface,  and  the  ambient. 


Fig.  2.  A  schematic  sketch  of  the  model  for  low  temperature  oxidation. 


In  the  low  temperature  regime,  the  flux  of  oxygen  through  the  external  layer,  must  equal  the  flux  of 
oxygen  through  the  porous  zirconia  channel.  Noting  that  oxygen  transport  occurs  as  permeation  of 
molecular  oxygen  through  liquid  boria  in  the  channels  that  is  driven  by  a  partial  pressure  gradient, 
the  flux  is  given  as 
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where  J  is  used  to  denote  flux,  f  is  the  fraction  of  pores  in  the  zirconia,  n  01_Bim  is  the  permeability 

of  oxygen  in  boria[20].  Rewriting  equation  (1),  we  obtain  the  following  equations  for  the  ratio  of 
the  layer  thicknesses  and  the  oxygen  partial  pressure  at  the  zirconia/extemal  boria  interface. 


T  (p^  -ps  ) 

^  _  j-  y1  02  1  02  J 


h 


(Pc 


02 


pZB 

‘  / 


02 


) 


or 


T  Pa 

pZB  _  ^r02 

r02 


+  fhatPo2 


ft1  ext  +L 


(2) 


cr 

m 

n 

7? 

□ 


Using  the  mass  balance  for  the  fonnation  rates  of  boria  and  zirconia,  and  the  inward  the  flux  of 
oxygen,  and  combining  with  equation  (1),  we  obtain  the  differential  equation  describing  zirconia 
layer  thickness  growth. 
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where  M  and  p  refer  to  molecular  weight  and  density.  The  rate  of  formation  of  boria  as  the  external 
scale  due  to  oxidation,  dnextBMj  /  dt  (=T i5203)  in  mole  per  unit  area  per  unit  time,  is  given  by  the  rate 

of  formation  of  boria  at  the  substrate  less  that  which  resides  in  the  porous  zirconia  and  that  which 
evaporates  from  the  surface. 
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At  low  temperatures  where  boria  is  too  viscous  to  flow,  the  thickness  evolution  of  the  external  boria 
scale  is  given  by  Equation  4(b).  At  temperatures  where  the  viscosity  of  the  boria  is  low,  significant 
flow  of  boria  will  occur  and  this  will  detennine  the  external  thickness,  as  given  by  the  following 
equation.  [21] 
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where  lspec  is  the  vertical  dimension  of  the  sample,  g  the  acceleration  due  to  gravity,  r/B203  viscosity 

of  boria  (data  from  [22]),  (j)  the  angle  the  face  subtends  to  the  horizontal  (90°  for  suspended 
specimens).  The  rate  of  evaporation  of  boria  at  the  surface,  under  boundary-layer-limited  diffusion 
is  given  by  the  following  set  of  equations  that  include  the  flow  velocity  of  the  ambient  gases  V,  its 
viscosity  r\,  gaseous  diffusion  coefficient  D,  and  specimen  length  Lspec. 
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Following  Graham  and  Davis  [23],  the  effect  of  the  boundary-layer-limited  diffusion  can  be  related 
to  the  maximum  evaporation  rate  given  by  the  Langmuir  equation,  through  a  parameter,  a,  to 
describe  the  evaporation  efficiency  under  the  actual  laboratory  condition. 
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Finally,  the  weight  change  per  unit  area,  AW/A,  and  recession,  R,  are  given  by: 
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The  differential  equations  (3)  and  (4)  are  solved  numerically,  and  used  to  predict  weight  change, 
recession,  and  external  scale  thicknesses. 

High  temperature  regime.  The  high  temperature  regime  has  a  much  simpler  structure  for  the  scale, 
as  shown  in  Fig.  3.  The  zirconia  thickness  is  taken  to  be  L,  the  superscripts  s  and  a  refer  to  the 
substrate  and  ambient.  The  partial  pressures  of  boria  and  oxygen  at  the  substrate/oxide  interface  are 
related  through  the  following  equation. 
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lZr02 


reaction®  ZrB2  J 
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,2/5 


(9) 


The  flux  of  oxygen  and  boria,  and  the  balance  of  fluxes  give  the  following  relationships. 
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Combining  equations  (9)  and  (10),  gives  the  following  polynomial  for  oxygen  concentration  at  the 
substrate,  CSQ  . 
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At  these  high  temperatures  (>1800C),  the  vapor  pressures  of  boria  are  so  high  thatC^,  =0  is  a  good 
assumption.  The  scale  growth  and  recession  are  then  given  by: 
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Fig.  3.  A  schematic  sketch  of  the  model  used  for  high  temperature  regime. 


The  evaporation  rate  of  zirconia  at  high  temperatures  cannot  be  neglected.  Using  thermodynamic 
data  on  vapor  pressure  of  ZrC>2  from  Barin  [24]  the  evaporation  rate  under  boundary-layer-limited 
diffusion  conditions  is  given  as: 


^/rO,  ~  ^ZrOn 


PZr0/(RT)-C 


Zr02 


=  a- 


1  ZrO, 


2  nMzroRT 


PZr0,atm  =  3x1 O10  cxp(--^^) 


(15) 


Comparison  of  Model  Predictions  with  Experimental  Data 

There  are  limited  data  in  the  literature  on  the  low-temperature  and  high-temperature  regimes  of 
oxidation  of  the  diborides.  For  the  low  temperature  data,  the  weight-gain  measurements  of  Tripp 
and  Graham  [8]  conducted  by  TGA  using  0.033  MPa  of  pure  oxygen  were  considered  the  most 
appropriate  data.  The  data  of  parabolic  rate  constant  by  Kuriakose  and  Margrave[25]  at  0.099  MPa 
Oo  was  also  considered.  For  the  high-temperature  oxidation,  the  only  data  available  are  the 
recession  data  summarized  by  Fenter[7]  from  the  work  of  Clougherty  et  al.[l,  2,  6].  Some 
calculated  model  parameters  that  result  from  the  experimental  conditions  of  quoted  data  are  shown 
listed  in  Table  1. 

Table.l  Experimental  conditions  and  resulting  model  parameters  at  1273K 

Diffusivity 


Reference 

Environment 

velocity 

sample  length 

Temp  (K) 

(B203) 

5 

a 

(m/s) 

(m) 

(K) 

(m2/s) 

(m) 

Tripp, Graham 

250Torr  02 

0.005 

0.02 

1173 

1.37E-04 

0.0424 

2.17E-05 

Kuriakose,  Margrave 

740Torr  02 

0.09 

0.005 

1173 

7.44E-05 

0.0034 

1 .47E-04 

Clougherty/Fenter 

air 

0.01 

0.01 

1173 

7.52E-05 

0.0144 

3.49E-05 

Berkowitz-Mattuck 

37.5T  02-He 

0.031 

0.003 

1173 

7.10E-05 

0.0044 

1.08E-04 

First,  we  show  the  model  prediction  across  the  entire  temperature  regime  to  illustrate  the 
temperatures  of  transition  and  how  the  transition  affects  several  parameters.  Fig.  4  shows  the 
predictions  for  zirconia  scale  thickness,  diboride  recession,  the  external  boria  thickness  and  internal 
boria  thickness  in  a  single  plot,  as  a  function  of  temperature  from  773K  to  2500K.  The  predictions 
are  for  the  conditions  used  by  Tripp  and  Graham,  i.e.,  0.033  MPa  of  pure  oxygen,  and  for  a  duration 
of  1  hour.  At  temperatures  below  800K  external  boria  forms,  but  above  800K,  the  flow  of  glassy 
boria  is  found  to  be  important.  The  model  predicts  that  this  external  boria  thickness  decreases  with 
increasing  temperature  and  is  non-existent  above  1273K.  Above  this  temperature,  the  boria  starts  to 
recede  inwards  and  is  all  gone  at  around  2073  K.  Thus  the  intermediate  temperature  regime  is 
predicted  to  be  between  1273K  and  2073K.  Note  that  the  recession  and  zirconia  scale  thicknesses  consult 
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are  not  affected  significantly  by  the  external  boria  layer  formation,  (<800K).  Obviously,  oxygen 
penneation  through  the  porous  channel  is  more  difficult  than  through  the  continuous  external  layer. 
Above  2098K,  there  is  a  brief  regime  where  Zr02  evaporation  is  negligible,  but  above  2400K,  it  is 
significant  and  accompanied  by  rapid  acceleration  of  recession. 


1/Temperature,  K 


Fig.  4.  Model  predictions  for  zirconia  scale  thickness,  diboride  recession,  the  external  boria 
thickness  and  internal  boria  thickness  are  shown  as  a  function  of  temperature  from  833K  to  2500K. 
The  predictions  are  for  an  exposure  time  of  1  hour  in  0.033  MPa  of  pure  oxygen. 


Fig.  5  (a)  shows  the  model  predictions  compared  with  the  data  of  Tripp  and  Graham[8]  for 
isothermal  weight-gain  versus  time  obtained  in  the  low  and  intennediate  temperature  regimes.  Fig. 
5(b)  shows  the  parabolic  rate  constants  measured  by  Tripp  and  Graham[8,  25]  at  0.033  MPa  of 
oxygen,  and  the  rate  constants  reported  by  Kuriakose  and  Margrave[25]  at  0.099  MPa  of  oxygen. 

There  is  reasonable  agreement  between  the  data  and  model  predictions.  Fig.  6  shows  a  comparison 
of  model  with  reported  measurements  of  recession  by  Clougherty  and  Fenter  [6,  7]  in  air,  and  by 
Berkowitz-Mattuck[5]  in  0.005  MPa  02-He  gas.  Good  correspondence  is  seen  except  for  the 
temperature  at  which  the  recession  accelerates  in  the  high  temperature  regime.  The  rate  of 
evaporation  of  zirconia  is  predicted  to  accelerate  only  above  2450K,  while  acceleration  is  seen  at 
around  2050K.  Spallation  of  the  scale  is  one  possible  explanation  for  this  discrepancy. 
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Fig.  5.  Comparison  of  model  predictions  with  experimental  data  on  (a)  isothermal  weight  gain  taken  from 
Tripp  and  Graham  [8]  (b)  parabolic  rate  constants,  taken  from  Tripp  and  Graham  [8],  Kuriakose  and 
Margrave  [25]  (model  shown  as  lines,  and  data  as  points) 


Fig.  6.  Comparison  of  model  predictions  with  experimental  data  on  recession  of  diboride.  The  high- 
temperature  data  in  air  (lh)  are  from  Fenter  [7],  and  the  intermediate  temperature  data  in  0.005  MPa  CF-FIe 
are  from  Berkowitz-Mattuck.[5]  (model  shown  as  lines,  data  as  points) 
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Summary 

The  transition  in  oxidation  mechanisms  of  diborides  was  modeled  as  a  function  of  temperature.  The 
model  assumed  that  at  all  temperatures  a  porous  zirconia  forms.  At  lower  temperatures  there  is  an 
external  layer  of  boria,  but  its  thickness  can  decrease  due  to  either  viscous  flow  or  evaporation.  In 
the  intermediate  temperature  regime,  the  boria  recedes  within  the  pores  of  zirconia.  At  the  high 
temperature  regime,  the  zirconia  is  dry  and  can  evaporate  as  gaseous  zirconia.  This  model  was 
evaluated  using  data  on  ZrB2  measured  under  different  experimental  conditions.  The  agreement 
with  experiments  is  found  to  be  very  reasonable  for  most  conditions  in  the  low  and  intermediate 
temperature  regimes.  The  temperature  at  which  zirconia  vaporization  accelerates  recession  in  the 
high  temperature  regime  is  overpredicted  by  several  hundred  degrees.  The  possibility  of  scale 
spallation  was  ignored  in  the  model,  but  may  have  been  present  under  the  experimental  conditions. 
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